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Lagrangian systems with nonholonomic constraints may be considered as singular differen- 
tial equations defined by some constraints and some multipliers. The geometry, solutions, 
symmetries and constants of motion of such equations are studied within the framework 
of linearly singular differential equations. Some examples are given; in particular the well- 
known singular lagrangian of the relativistic particle, which with the nonholonomic constraint 
v 2 = c 2 yields a regular system. 
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X. 

c3 ■ 1 Introduction 

The main goal of this paper is to study the nonholonomic mechanical systems within the frame- 
work of linearly singular differential equations. 

Nonholonomic mechanical systems, i.e., mechanical systems with non-integrable kinematic 
constraints, have been discussed since the last years of nineteenth century. However, the geomet- 
ric foundations for the theory were given in |VF 72j . Since then, several approaches have been 



taken to deal with the subject, for instance, a hamiltonian approach in BS 93 , a lagrangian 
approach in LM96j, a more general Poisson framework in |Mar 98j . or an approach based on a 
gauge independent formulation of lagrangian and hamiltonian mechanics in |MVB 02] . Symme- 
tries of these systems, as well as reduction schemes derived from them, have also been considered 
in the literature, see jKoi 921 IHKMM 961 IK M 981 ICL991 1 Mar 03] . 

A lagrangian system with nonholonomic constraints may be considered, more generally, as 
a singular differential equation defined by some constraints and some multipliers: 

x = g(x) + ^■i/ t /i M (x), (f) a (x) = 0. 
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Such an equation can be described geometrically as a linearly singular differential equation, 
that is, a differential equation where the velocities are not isolated because of a linear factor 
multiplying them: 

A(x)x = b(x). 

This is a special type of implicit differential equation. The idea of modelling mechanical systems 
as implicit differential equations is found in earlier papers by Tulczyjew |MT 781 IMMT95] . and 
it has also been used to deal with nonholonomic constraints |Tul 86| [ILMM 96 . 

Linearly singular differential equations were geometrically presented in |GP 91| IGP 92j . This 
general framework includes for instance the presymplectic systems and the lagrangian formalism 
GP 92] , the higher order singular lagrangians and their "higher order differential equation" 
conditions |GPR9lj . as well as many other systems that appear in technological applications. 
To solve the corresponding equation of motion a consistency algorithm can be performed. This 
algorithm is indeed a generalization of the presymplectic constraint algorithm |GNH 78] . 

We will see that a system with constraints and multipliers, and in particular any nonholo- 
nomic mechanical system, can be described as a linearly singular system. This implies that all 
the methods and results about these systems can be applied directly to nonholonomic systems. 

More precisely, the combination of two operations that can be performed on linearly sin- 
gular systems — restriction to a subsystem and projection to a quotient — can be applied to 
obtain what we call a generalized nonholonomic system. In particular, we discuss the regularity, 
consistency and equations of motion of these derived systems. 

The symmetries of a linearly singular differential equation have been studied in |GP02j . In 
this paper we consider the relation between the symmetries of a system with nonholonomic 
constraints and the symmetries of its original unconstrained system, both modelled on linearly 
singular differential equations. We also study their constants of motion. 

The paper is organized as follows. In section 2 we give some definitions and results regarding 
linearly singular differential equations, their solutions and their symmetries. In section 3 we 
introduce generalized nonholonomic systems and discuss some of their properties. Symmetries 
and constants of motion of generalized nonholonomic systems are discussed in section 4. In 
section 5 we show how a lagrangian system with nonholonomic constraints can be described in 
terms of a generalized nonholonomic systems. The case of a relativistic particle is studied in 
section 6, where we see that a nonholonomic constraint can convert a singular lagrangian into a 
regular system. Two additional examples are studied in section 7. Finally, an appendix contains 
some auxiliary results formulated within the framework of linear algebra. 

2 Previous results: linearly singular systems 

In this section we recall some definitions and results from |GP 911 IGP 921 IGP 02| . 

Let M be a manifold. An implicit differential equation on M is defined by a submanifold 
D C TM. A path £: / — ► M is a solution of this equation when 



ZV) C D. 
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In coordinates, if the submanifold D is described by some equations F = and the path £ is 
represented by some functions x(t), then the local expression of the implicit differential equation 
is F(x, x) = 0. 

We have a particular case when D = X(M), with X a vector field on M. Then X defines 
an explicit differential equation, and £ is a solution iff 

e = Xo£. (2.2) 

Now the local expression is x = f(x). 

A linearly singular differential equation on M is defined by a vector bundle it: F — ► M, a 
vector bundle morphism A: TM — > F, and a section f\M —* F of 7r. A path £: I — > M is a 
solution when 

Ao£ = /o£, (2.3) 

whose local expression is ^4(x)x = /(x), with A(x) a (singular, in general) matrix. 

We denote by (A: TM — > .F, /) the linearly singular system. The following diagram shows 
all these data: A > ^ 





- M 

The associated implicit differential equation is 

D = A~ 1 (f(M)) C TM. (2.4) 

We say that the linearly singular differential equation is regular when A is a vector bundle 
isomorphism. In this case, the associated explicit differential equation is given by the vector 
field X = A~ 1 o f. 

The solutions of the system can be equivalently described as integral curves of vector fields. 
Let us remark that in general the solutions are restricted to a submanifold S C M because 
the equation Q2.3|) may not have solutions passing through every point x 6 M. Therefore, the 
equation of motion can be written as an equation for a vector field X and a submanifold S: 

{X tangent to S 

where the notation means equality at the points of S. 

A recursive algorithm can be applied to find the solutions of a linearly singular differential 
equation. Its first step is to note that, in order that a solution passes through a point x G M, 
it is necessary that 

f(x)elmA x , (2.6) 
so the solutions are necessarily contained in the primary constraint subset 

Mi = {x G M | f{x) G ImAJ, (2.7) 
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which will be assumed to be a closed submanifold. The tangency to M\ forces the initial system 
to be restricted to (Af.TMi — > where Ay = -A| TM , F\ = F\ M and fy = f\ M . The 

algorithm follows recursively, and, under some regularity assumptions at each step, it ends with 
a final constraint submanifold 5* such that f(S) C Im As; thus the system is consistent, and the 
equation 

A s oX = f s (2.8) 

for a vector field X tangent to S has solutions. Given a particular solution X , the set of 
solutions of (|2.5|) is X + Ker As- 

We finish this section by giving some definitions and results about symmetries. A symmetry 
of a linearly singular system (A: TM — > F, f) is a vector bundle automorphism (ip, <I>) of tt: F — > 
M such that 

/ = $*[/]:= $o/o yT 1 , A = $4i]:=$oio(T^)- 1 . (2.9) 

An infinitesimal symmetry of a linearly singular system (A: TM — > i 7 , /) is an infinitesimal 
automorphism (V, W) of the vector bundle ir: F — > M such that its flow (Fy , F|^) is constituted 
by local symmetries of the linearly singular differential equation. The last property is equivalent 
to the following conditions: 

TfoV = Wof, TA o V T = W o A, (2.10) 
which are the infinitesimal version of (|2.9|) . 



3 Generalized nonholonomic systems 
The geometric setting 

Among the various operations that can be performed with a linearly singular system {B: TN — > 
G,g), we are especially interested in the subsystem defined on a submanifold j:M N, and 
the projection p:G—> G/G' to a quotient with respect to a vector subbundle G' C G: 




N M N 



Suppose that the original system admits solutions Y on a submanifold Nf C N. Then the 
subsystem on M has solutions on the submanifolds of M D A r f over which a solution K of the 
initial system is tangent. On the other hand, the quotient system has, in general, more solutions 
than the initial system: if Z is any vector field on tangent to JVf with values in i? -1 (G") then 
Y + Z is a solution of the quotient system on Nf, there may also exist solutions defined on a 
submanifold greater than Nf. 

It is well known that the dynamics of systems with nonholonomic constraints is a mixture 
of both constructions: the presence of some constraints, combined with a certain degree of 
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arbitrariness expressed through some multipliers. This combination may result advantageous: 
though in general Y is not tangent to the submanifold M, it may happen that for some vector 
fields Z in B~ 1 (G') one has solutions Y + Z tangent to M, or at least to a "big" submanifold 
of M. 

In this paper we will call a generalized nonholonomic system the linearly singular system 
(A:TM — > F, /) defined from (B:TN — > G, g) by a constraint submanifold M C N and a 
subbundle of constraint forces G' C G| M as follows: 

. F = (G| M )/G', 

o 

• A=po B| M 0T7, and 

• f =P° g\ M > 

o 

where p: G| M — > is the projection to the quotient, and Tj denotes the tangent map 

of j with the image restricted to M. All this is shown in the following diagram: 



A 



TM 




M = M 



Regularity and consistency 

Before discussing the equations of motion, we want to study some general properties of the 
generalized nonholonomic system (A: TM — ► F, /), namely, whether A is surjective (we will also 
say that the system is surjective) or bijective (the system is regular), or the equation AoX = f 
is everywhere consistent. 
Let us denote 

H = B- 1 (G')CTN\ M , 
which is a vector subbundle whenever the morphism B has constant rank. 

Proposition 1 With the preceding notations, the generalized nonholonomic system is surjective 
iff 

B(TM) + G'= G\ M . 

Assuming that the original system is surjective, the nonholonomic system is surjective iff 

TM + H= TN\ w, 



and it is regular iff in addition 

TMHH = {0}. 
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Proof. We want to decide whether A = p o B\ M o Tj (the composition of an inclusion, a 
morphism and a projection) is surjective or injective, and this is given by lemma 1 in the 
appendix. ■ 

The preceding result could be refined also in the case where B is injective, but this does not 
seem so interesting. As an immediate consequence, we have: 

Corollary 1 Suppose that the original system is surjective (or, more particularly, regular). 
Then the generalized nonholonomic system is regular iff 

TN\ M = TM © H. 



These relations can be given a more concrete form in terms of constraints and frames. 
Consider a local basis n)i<ft<m of sections for the subbundle H C TN\ M (they are vector 
fields in N, but defined only on M). Consider also a set of a Q constraints <p a , linearly independent 
at each point, that locally define the submanifold M C N. Finally, consider the matrix 

D- = (d^| M ,r /i ) = r /i -^, (3.i) 

whose elements are functions on M. 
Proposition 2 With the preceding notations, 

1. TM n H = iff rank(L>^) = m . 

2. TM + H= TN\ M iff rank(D^) = a . 

3. TM © H = TN\ M iff (D°£) is a square invertible matrix. 

Proof. It is a consequence of lemma 3 in the appendix, since the d(j) a \ M constitute a basis for 
the annihilator of TM in (TiV| M )*. ■ 

The connection of such a matrix with the notion of regularity and consistency of a constrained 
system was already noted in |CR93| ILM 96j . 



Equations of motion 

From the definition of the generalized nonholonomic system (A:TM — > F, /), it is clear that a 
path £: / — > TV" is a solution of the equation of motion iff it is contained in M and 

fio(- 9 o(eG', (3.2) 
If some sections A u constitute a frame for G' , then this equation can be written as 

Bot = got + Y,v v &v°€, (3-3) 

for some multipliers v u (t). 
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In the same way, for a submanifold S C M and a vector field X on M tangent to S, the 
equation of motion iol^ / can be written as 

BoX-geG', (3.4) 

where the equation must only hold on the points of S. This equation may be also written as 

BoI-g + ^/A„, (3.5) 

s 

v 

for some multipliers v u {x). 

Of course, we can apply the constraint algorithm to find the solutions of this linearly singular 
system. However, there is an alternative way to solve the problem when the original problem is 
regular, or at least consistent. Under this hypothesis, let Y be a vector field on N, solution of 
the equation of motion of the linearly singular system {B: TN — > G,g): 

B o Y = g. 

(For most applications the original system is regular, and then the unique solution of this 
equation is the vector field Y = B^ 1 o g.) 
Using Y, the equations of motion become 

i-Yo^eH (3.6) 

for a path £ in M, and 

X - Y C H, (3.7) 

s 

for a vector field X on M that should be tangent to S. 

These equations can be expressed in a more concrete form in terms of the local basis (r^) 
of sections for the subbundle H C TN\ M : 

e = yo£ + ^uT M o£, (3.8) 



for some functions u^(t), and 



X^Y + ^uTp, (3.9) 



for some functions on M. 

Let us examine whether this last equation has solutions. The requirement for X of being 
tangent to M is X ■ (fr a ~ 0, which reads 

M 

J2D a ^ + Y-r^0, (3.10) 

where (D") is the matrix defined by ()3.1|) . From this it is clear that the generalized nonholonomic 
system is regular iff the matrix (-D") is invertible on M, and in this case the equation H3.1Uj) 
directly determines the functions that give the solution X expressed in (|3.9fl . More generally, 
the nonholonomic system has solutions if the matrix (D") has rank a a . 

Geometrically, the decomposition TN\ M = TM®H stated in Corollarynhas two associated 
projectors V, Q- Writing Y = V oY + QoY on M, the following result is clear: 
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Proposition 3 With the preceding notations, if the original system is consistent, with a solu- 
tion Y , and the generalized nonholonomic system is regular, with solution X, the latter can be 
obtained as 

X = VoY\ M . (3.11) 

■ 

Such projectors were studied, in the context of nonholonomic lagrangian systems, in [LM 96 . 

4 Symmetries and constants of motion 

Let us consider a generalized nonholonomic system (A:TM — > F,f), obtained from a linearly 
singular system (B:TN — > G,g) by means of a restriction to a submanifold M C N and a 
projection to the quotient p: G\ M — > (G\ M )/G' , where G' C G\ M is a vector subbundle. 

Recall the definitions of symmetry and infinitesimal symmetry given in section 2. Our aim 
is to study the relation between the symmetries of the original linearly singular system on N 
and the symmetries of the generalized nonholonomic system on M. In the next proposition, we 
give sufficient conditions on a symmetry of the original system in order to define a symmetry of 
the constrained system: 

Proposition 4 Let be a symmetry of (B:TN — > G,g). Suppose that ip leaves the sub- 

manifold M C N invariant, and ^ leaves the subbundle G' C G\ M invariant. Then (99,$), 
where ip = ip\ M , and <3?: (G\ M )/G' — > (G\ M )/G' is the map induced on the quotient from \P, is 
a symmetry of (A: TM — > F, /). 

Proof. We have 

A o Tip = p o B oTj o T( V>Im) = P° B o Tip o Tj = p o iff o B o Tj = $ o p o B o Tj = $ o A, 
and 

foip=pogo i/j\ M =po$ o g\ M = $opo g\ M = $ o /, 

so the two conditions for being a symmetry are satisfied. ■ 
We can obtain a similar result for infinitesimal symmetries, by making use of their infinites- 
imal characterization ()2.1Uj) : 

Proposition 5 Let (V,V) be an infinitesimal symmetry of (B:TN — > G,g). Suppose that V is 
tangent to the submanifold M C N, and V is tangent to the subbundle G' C G\ M . Then (U, U), 
where U = V\ M and U:(G\ M )/G' — > T{{G\ M )/G') is the vector field induced on the quotient 
from V , is an infinitesimal symmetry of (A: TM — > F,f). 

Proof. The proof runs as in proposition 21 

Tf o U = Tp o Tg o V\ M = TpoV o g\ M = U o p o g\ M = U of, 

TAoU T = Tp o TB o T(Tj) o | TM = Tp o TB o V T o Tj = 
= TpoyoBoTj = [/opoBoTj = {/oA 
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We now consider constants of motion. Suppose that the original system has a solution 
Y G X(N), and let us consider a function h G C°°(N) such that Y ■ h = 0. Under which 
conditions is h\ M a constant of motion of the generalized nonholonomic system? 

Suppose that both the original system and the nonholonomic system are regular, so that 
TN\ M = TM © H; let V be the projector to the first factor, which, according to Proposition EJ 
relates the dynamics of both systems as X = V o Y , Then we have a simple characterization: 

Proposition 6 With the preceding hypothesis, write X = Y — V , where V is a section of H C 
TN\ M . Let h be a constant of motion of the unconstrained system. Then h\ M is a constant of 
motion of the generalized nonholonomic system iffT ■ h = 0. 

Proof. It is straightforward: 

X - h = (Y - V) -h = Y -h-Y -h. 
(Note that Y and T, considered as sections of TN\ M , map functions on N to functions on M.)m 

5 Lagrangian systems with nonholonomic constraints 

In this section we will show that the dynamics of a lagrangian system with nonholonomic con- 
straints (the nonholonomic mechanics) falls into the class of generalized nonholonomic systems 
of section 3. 

We begin by considering a configuration manifold Q, its tangent bundle TQ, and a lagrangian 
function L: TQ — > R. The lagrangian mechanics may be described as the linearly singular system 
(0 L :T(TQ)^T*(TQ),d£ L ). 

T(TQ) — — ^— T*(TQ) 




TQ 

Here El is the lagrangian energy and ul is the Lagrange's 2-form. Though we do not want to 
dwell on these well known objects, some properties of lol and the vertical endomorphism will 
be needed later, so let us briefly recall them. See |Car 90j for more details. 

First, we have the vertical endomorphism J of T(TQ), whose kernel and image are the vertical 
subbundle V(TQ). Its transposed morphism is an endomorphism J of T*(TQ), whose kernel 
and image are Sb(TQ), the bundle of semibasic forms. This is used to define the Lagrange's 
1-form Ol = * J o dL and 2-form u>l = —d0L on TQ. 

From now on we consider the case of the lagrangian being regular, which amounts to ujl 
being a symplectic form. Then, it induces a vector bundle isomorphism u)x,:T(TQ) — > T*(TQ) 
mapping vertical vectors to semibasic forms, thus yielding an isomorphism V(TQ) Sb(TQ). 

It is well known that the lagrangian dynamics on TQ is described by the only vector field 
Xl solution of 

Q L ° X L = dE L - 
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Note that it is a second-order vector field. 

j 

Now let us introduce the nonholonomic constraints, which define a submanifold M TQ of 
dimension m. We will consider only the case where this submanifold restricts the velocities, not 
the configuration coordinates. In a more formal way, this is described by the conditions given 
in the next proposition: 

Proposition 7 Let M C TQ be a submanifold. The following conditions are equivalent: 

1. The projection M — > Q (restriction of the tangent bundle projection tq:TQ — > Q) is a 
submersion. 

2. (TM) h nSb(TQ)| M = 0. 

3. The submanifold M C TQ can be locally described by the vanishing of some constraints (j) 1 
whose fibre derivatives J-(j) 1 are linearly independent at each point of M . 

4- The submanifold M C TQ can be locally described by the vanishing of some constraints (j) 1 
such that the 1-forms A* = *J o d(f> 1 are linearly independent at each point of M. 

( dcf? \ 

In coordinates, these conditions mean that — -=- has maximal rank. ■ 

\dv k J 

Note that under the preceding conditions the image tq{M) C Q is an open submanifold, 
and so we can replace Q with this submanifold. So, from now on, we assume that the projection 
M — > Q is a surjective submersion. 

Now we will consider the following vector bundles: 

TM C T(TQ)| M , 

{TMf C T*(TQ)| M , 
G> := *J((TM) h ) C Sb(TQ)| M , 
H := u-\G') C V(TQ)| M . 

Suppose that M C TQ is defined by the vanishing of some independent constraints (ff as in the 
preceding proposition. Then (TM) h is spanned by the d(/> l | M . We denote by A 1 and T{ their 
corresponding images in G' (through *J) and H (through S^ 1 ). The following diagram shows 
all these objects: 

TMw T(TQ)\ M — T*(TQ)|^ {TMf = (d^| 



\M> 



J 



(T i )=H^ V(TQ)| M -^L^ Sb(TQ)| M ^ G> = (A 4 ) 

So we have two subbundles TM, H C T(TQ)\ M . We have rankTM = m and rank(TM) h 
n — m; the conditions in Proposition [7| imply also that r&nkH = rankG' = n — m. 
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Theorem 1 The nonholonomic mechanics defined by the lagrangian L and the constraint sub- 
manifold M C TQ is the generalized nonholonomic system defined from the lagrangian mechan- 
ics (Ql:T(TQ) — > T*(TQ), d£x) by the constraint submanifold M C TQ and the subbundle of 
constraint forces G' = *J((TA/) h ) C T*(TQ)\ M . 

o 

TM T(TQ)| M _^U T *(TQ)| M * T*(TQ)| M /G' 




M — M 
Proof. The equation of motion for a path £ = 7 such that G M is 

e = ^o£ + £yr i0 £. (5.1) 

Instead, let us write the equations of motion for vector fields: according to (|3.5j) . for a second- 
order vector field X on TQ, tangent to M, the equation is 

ijcWi^dEx + XVAi. (5.2) 

i 

or, according to (|3.9|) . 



X^X L + ^nTi. (5.3) 



But in coordinates equation (|5.2j) reads 



d /dL\ _ ^ i< 
9g dt \dv J ^-r 1 dv 



which is the equation of motion of the nonholonomic mechanics defined from L and the con- 
straints — see for instance |Arn 83j . ■ 

If, in addition to (TM) h n Sb(TQ)| M = 0, we have TMnH = 0, then T(TQ)\ M = TM®H, 
and so there is a unique solution X of the equation of motion, which can be obtained from Y 
through the projector to TM as described by Proposition OJ 

The case of a singular lagrangian 

The preceding method can be conveniently adapted if the lagrangian is singular. Of course, 
one can not use the direct sum decomposition. However, the formulation of the nonholonomic 
dynamics as a quotient system on a submanifold remains unchanged, except that the second- 
order condition is not automatically satisfied by X and must be imposed as an additional 
equation for it: 

J o X ~ Ato- 

m y 

This condition may be included in the equation of motion of the nonholonomic dynamics in 
the same way as can be done with the lagrangian dynamics, using the time-evolution operator 
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K of lagrangian dynamics BGPR86 |GP 89) . With it, the lagrangian dynamics is the linearly 
singular system „ 

T(TQ)^S)TQx^ L T(T*Q) 




TQ 

where TL:TQ — ► T*Q is the Legendre's transformation (fibre derivative) of L. In terms of 
vector fields, the lagrangian dynamics is thus defined by the equation 

T{TL) o X ~ K. 

Then, it is readily seen that the nonholonomic equation of motion can be written 

T(fL)oI~K-^u¥'. (5.4) 

Here is a certain vector field along J-L, which is defined from the fibre derivative of a function 
</>: T*Q -> R —see |<;P01| for details. 



6 Relativistic particle with a nonholonomic constraint 

In this section we study the motion of a relativistic particle as a nonholonomic constrained 
system. We will consider two possible lagrangian functions, a regular one (deeply studied 
in |KM01| ) and a singular one. 

Let us consider a particle with mass m and charge e moving in spacetime. We model 
spacetime as a 4-dimensional manifold Q, endowed with a metric tensor g of signature (1,3). 
Suppose furthermore that the particle is subject to the action of an electromagnetic field F = dA, 
where A G O x (Q), and a potential U G C°°(Q). 

Recall that there are some relevant objects associated with the metric g, namely, the isomor- 
phism g: TQ — > T*Q (we will denote X' = goX)), the Levi-Civita connection V, the differential 
forms 6 g = g*{9 Q ) G ^(TQ) and u g = g*{uJ Q ) = -dO g G ^ 2 (TQ), the energy E g {u q ) = 
^g(u g ,u q ) G C°°(TQ), and the geodesic vector field S g , which satisfies is g u g = dE g . We denote 
v = ^2E~. 

We will study two different lagrangian functions, namely 

Li(u q ) = -mcg(u q ,u q ) 1/2 - ^(A(q),u q ) - U(q), 

and 

1 e 
L 2 (u q ) = --mg(u q ,u q ) - -(A(q),u q ) - U(q). 
Zj c 

Forgetting the potential, L\ is the singular lagrangian commonly used in relativistic mechanics 
to describe a particle in an electromagnetic field; it is defined only on the open set of time-like 
vectors of TQ. The lagrangian L 2 appears in |KM01| . Our aim is to compare both systems, 
and to introduce the nonholonomic constraint v 2 = c 2 to them. 



X. Gratia and R. Martin, "Regularity and symmetries of nonholonomic systems" 13 



mc i 



The lagrangians L\ and L 2 have, respectively, associated Lagrange's 1-forms 6\ — — —w g - 
- c TqA and 6> 2 = -m9 g - \tqA; the Lagrange's 2-forms are u)\ = -^uJg - ^dv A9 g + \tqF and 
u>2 = —mujg + ^TqF; and the lagrangian energies are E\ = U and E2 = — \mv 1 + U. 

The symplectic formulation of the equations of motion for the lagrangians L\ and L2 are, 
respectively, 

i X 0Ji = dEi, (6.1) 

and 

i x uj 2 = dE 2 , (6.2) 

for second-order vector fields X. For any 2-form cj, we will also denote ixu by uj{X). 

It is worth writing down the Euler-Lagrange equations of motion for a path 7, which are, 
for lagrangians L\ and L 2 : 

and 

"i(V t 7) b + -UfF - dU = 0. (6.4) 
Let us now consider equations (|6.1|) and (|6.2j) . 

As uji is not surjective, equation (jfi.lj) could have no solutions. We denote by A = (f ^ the 
Liouville vector field, T = q l J=^ the natural vector field along tq, and £ v the vertical lift of a 
vector field £: TQ — > TQ along tq. We have that Keri^i = (A, £), where 

^ = S g -— 2 ({i T Ffy. (6.5) 
mr 

We can see that a)i(^(grad C/) v ) = d£7 — (^jiTdU)9 g and that # 9 ImcDi. Therefore 
equation H6.1|) has solutions if and only if i^dU = 0, that is, the potential U is constant, which, 
in practice, is the same as taking U equal to 0. 

Since S is a second-order vector field, in absence of potential the solutions of equation (|6.1|) 
are X\ = T, + fiA, where fi is an arbitrary function. If, in addition, there is no electromagnetic 
field, then the solutions are S g + /j,A, and their integral curves are reparametrized geodesies. 

On the other hand, equation (|6,2j) is regular, and its solution is 

X 2 = S g + -(grad[/) v - —{(i T Ff)\ (6.6) 
m mc 

This can be proved making use of the relations iz^oj g = —Tq(Z^) for vector fields Z along tq, 
and is(tqF) = Tq{itF). In this case, in absence of electromagnetic field and potential, the 
solutions are the geodesies of g. 

Now we introduce the nonholonomic constraint 

0(«<j) := g{u q ,u q ) - c 2 = 0, (6.7) 



which defines a submanifold M C TQ. 
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The subbundle of constraint forces is ( t J(d(j)))\ M = (0 g )\ M , therefore, according to equa- 
tion (|5.2|) . the equations of motion for both lagrangians become 



and 



ixoJi ~ d^i + X9„, (6.8) 

M 



i x uJ2 ^dE 2 + X9 g , (6.9) 

M 



for second-order vector fields X tangent to M. 

Note that if a path 7 satisfies the constraint then it also satisfies the equation = ~^9iii t) = 
25(7, Vt7), so looking at equations (|6.3j) and Q6.4|) we realize that the two constrained systems 
have the same equations of motion: 

(V t7 ) b H 
5(7. 7) ~ "~ 



m(V t7 ) b + -iif - dU = Xj\ 

c 7 ' ' (6.10) 



The multiplier A can be found by contracting the equation with 7, which gives A = — ^i^dU. 

We are going to see this equivalence of the solutions of both Euler-Lagrange equations by 
computing the solutions of equations (|6.8|) and (|6.9|) . 

First let us analyse equation Q6.9[) . From A • <j) = 2v 2 ~ 2c 2 7^ and ia^>2 = m @g, it follows 
that TM ( (0 g )\ M ) = (TQ)\ M , so, by proposition^ the system is regular. Its solution is 
X = X2 + ^A, where the multiplier A is found by imposing that X is tangent to M: 

O = X- = X 2 -0 + —A • 6 ~ — i T dU + 2— c 2 . (6.11) 

m Mm m 

Therefore, the solution of the second system is 



X = S g + -(grad U) v - —{(i T F)*) v - ^-(i T dU)A. (6.12) 



Now let us analyse equation l|6.8|) . Since y = ^(gradC/) v — -^-^(irdU)A is a vector field 
tangent to M and ^i(y) dJ7 — {-^iTdU)9 g , the system is consistent. We can see that 

TM n wf 1 ( (0 g )\ M ) = TM n Ker Qi = (£}| M , (6.13) 

so the system is not regular. Then, the solutions of the equation are Y + riY,. Since Y is vertical, 
in order to be a second-order vector field the function fj, must be equal to one, so the solution is 
Y + £ ~ X, exactly the same as for the lagrangian L 2 . 

M 

7 Examples 
Example 1 

Consider the differential equation on N = R 2 defined by the vector field Y = J| + yJ^. We 
restrict this system to a generalized nonholonomic one by means of the construction of section 3, 
taking the submanifold M = R X {a} C N and the subbundle C = (xJ% + J^) C TiV| M . 
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In this case TN\ 



M 



TM 



C and the projectors associated with this decomposition are 
V: 



d_ 

d.r 

a_ 

dy 



d_ 

d.r 



Q: 



8_ 

d.r 

a_ 

dy 



a_ 



x dx ' dy' 



Thus X = VoY 



M 



(1 



ax. 



JL\ 

dx I M 



is the solution of the generalized nonholonomic system. 



Let us study the infinitesimal symmetries of both systems. We can see that a vector 
field V € X{N) is an infinitesimal symmetry of the unconstrained system if it has the form 
V = V 1 (ye~ x )^ + eX V 2 {y e ~ x )~§[;i where V 1 and V 2 are arbitrary smooth functions. 

On the other hand, since the constrained system is one-dimensional, its infinitesimal sym- 
metries are the vector fields U = kX, with k £ R. Observe that, in principle, an infinites- 
imal symmetry of Y does not lead to an infinitesimal symmetry of X by restriction to M, 
even when Y\ M G X(M). Nevertheless, if we also require that V T (C) C TC, then we obtain 
V l {t) = k(l + aln(t/a)) and V 2 (t) = 0, so that actually V\ M = k(l — ax) -§^\ M is an infinitesimal 
symmetry of X. 



Example 2 

Here we discuss an example of a particle with a nonholonomic constraint, due to Rosenberg 
|Ros 77j . Consider a particle moving in R 3 with lagrangian function 

1 



L = ^(x 2 + y 2 + z z ) 



subject to the nonholonomic constraint 



z — yx. 



Using the notation of section 5, we have N = TR 3 , u>l = dx A dx + dy A dy + dz A di and 
dE^ = xdx+ydy+zdz, so the unconstrained dynamics is the well-known free dynamics described 
by the vector field 

^ = 0- L \dE L ) = x§- x+ y^ + z§-. 
The constraint submanifold is M = {z=yx}, with tangent bundle 

TM 

and the vector subbundle C C T./V| M is 



[ ■ I / 9 d . d d d d d 

dx ' dy X dz ' dz ' dx ^ dz" 1 dy 



M 



c = {air . m)))=(v --- 



d 



d 



M 



Note that TN\ M splits as TN\ M = TM©C, so the only solution X of the constrained lagrangian 
system is the projection of Xj\ M to TM according to this decomposition: 



X 



. d . d . d 

dx ^ dy dz y 2 + 1 dx ' y 2 + 1 dz 



yyx d yx d 
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We choose (x,y, z,x,y) as coordinates on M. With this system, the vector field X reads as 

. d d . d yyx d 
X = x-t — h y— — h yx- 



dx dy dz y 2 + 1 dx 

After some calculus, we can find the symmetries and constants of motion of both systems. The 
constants of motion of the free particle are the functions G(x, y, z, xy — yx, yz — zy), where G is 
an arbitrary function with five variables. The infinitesimal symmetries are linear combinations 

r\ f\ f\ r\ Pi Pi Pi Pi Pi 

of the six vector fields ^, x-q^ + x m , y-^ + y^ and z-g^ + z m , with the constants of 
motion as coefficients. 

The constants of motion of the constrained system, written in coordinates of M, are 



F [y, x\Jy 2 + 1, yx - arcsmh(y)x^/y 2 + 1, yz - x(y 2 + 1) ) , (7.1 



and the infinitesimal symmetries are linear combinations of the five vector fields 

d d d d d xyy d 
— , — ,x— + y— + yx- 



dx ' dz ' dx dy dz y 2 + 1 dx ' 



argsinh(y) d \Jy 2 + 1 d ^ 1 d 



y dx y dz ^Jyi _|_ i dx ' 

(y — argsinh(y)y / y 2 + 1) d yd x d xy 2 d d 



y 2 dx y dy y 2 dz y(y 2 + 1) dx dy 

with the constants of motion as coefficients. 

In order to illustrate proposition we take a function g = G(x, y, z, xy — yx, yz — zy), i.e., a 
constant of motion of Xl, such that Z ■ g = 0, where Z is the section of C 



xy ( d d 
Z = X L\ M -X=— x \y--- 



M 



This yields to 

g = H (y, \J z 2 + x 2 , z + yx — xy — arg sinh(i / x) \J z 2 + x 2 , yz — zy — x^j 
and we see that g\ M is just the expression l|7.1j) . 



Appendix: some lemmas about linear algebra 

Here we collect some results about linear algebra on vector bundles that are needed in section 3. 
These lemmas are stated and proved for vector spaces, but of course nothing changes essentially 
if vector bundles are considered instead. 
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Lemma 1 Let f:E^>Fbea linear map between vector spaces, and E Q C E and F C F vector 
subspaces. Denote j: E — > E the inclusion, p:F—> F/F a the projection to the quotient, and 
consider the composition f = p o f o j . Then: 

1. f is injective iff E Q n / _1 (-F ) = {0}. 

Assuming f injective, this also amounts to f(E ) n F Q = {0}. 

2. f is surjective iff f(E ) + F = F. 

Assuming f surjective, this also amounts to E + f^ 1 (F ) = E. 

3. When f is surjective, f is bijective iff E a © f~ 1 (F ) = E. 
When f is injective, f is bijective iff f(E Q ) © F Q = F. 

Proof. First note that 

Kerf = E nf-\F ), Im / = (f(E Q ) + F )/F Q . (A.l) 

These equalities are clear: the kernel is constituted by the vectors in E Q mapped to F by /, 
and the image of a subspace F' C F by p is (F' + F )/F . This readily yields the first assertions 
about injectivity and surjectivity. 

Their equivalent formulations when / is injective [or surjective] can be proved using the 
formulas for f(E 1 D E 2 ) and f~ 1 {F 1 n F 2 ) [or for the sum], as well as f~ l {f{E )) = E + Ker /, 
f(f- 1 (F ))=F nlmf. 

Finally, the assertions about the bijectivity of / are a trivial consequence of the other ones.H 

Remember that a linear equation f(x) = b is consistent iff b E Im/. Now let us study a 
linear equation on E defined as in the preceding lemma by / and the class of an element b £ F. 

Lemma 2 The linear equation f(x) = b is equivalent to the couple of equations f(x) — b £ F Q , 
x € E a . It is consistent iffb G f(E ) + F ; in this case the solution is unique iff E Q n f^ 1 (F a ) = 
{0}. . 

Finally, let E C G be a subspace of a vector space. Recall that the annihilator (or orthogonal) 
of E is the subspace 

E h = { 7 <E G* [ (Vz G E) ( 7 , x) = 0} C G* . 

This space has a close relationship with G/E. Indeed, the transpose map of G — > G/E defines 
a canonical isomorphism 

5: (G/E)* ^E h , 
such that, for a G E h and z G G, (S^ 1 (a), z + E) = (a, z). 

Lemma 3 Let E,F C G be vector subspaces. Let (a 1 , . . . , dP) be a basis for the annihilator 
E^ C G* , and (v\, . . . ,v q ) a basis for F. Consider the matrix D = (Z? l )i< l < P with elements 

D) = (a % ,Vj). Then: 

1. E + F = G iffi&nkD=p. 
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2. E n F = {0} iffvankD = q. 

3. E © F = G iff D is square invertible. 

Proof. Consider the linear map e: F — > G/E denned as the composition of the inclusion F <^-> G 
and the projection to the quotient G -» G/E. It is clear that .E + -F = G iff e is surjective, and 
EPiF = {0} iff e is injective, so the only thing to prove is that the given matrix is the matrix D 
of e in appropriate bases: the basis (vj) for F, and the basis (cvj), the dual basis of a 1 = 5^ 1 (a t ), 
for G/E. 

Then, if e{yj) = oiiD 1 -, we have D l - = (a l ,£(vj)) = (a l ,Vj), which is what we wanted to 
prove. ■ 
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